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ABSTRACT 

We  present  a  low  resolution  (R  =  X  /  AA.  ~  200)  far- 
infrared  (43-197  pm)  spectrum  of  the  ultraluminous 
galaxy  Arp  220,  obtained  from  deep  full-range  scans 
using  the  Long  Wavelength  Spectrometer  (LWS)  aboard 
the  Infrared  Space  Observatory  (ISO).  The  spectrum  is 
vastly  different  from  the  spectra  of  less  luminous  IR 
bright  galaxies:  the  fine-structure  lines  that  are  typical 
of  IR  bright  galaxies  are  weak  or  absent  and  the 
spectrum  is  dominated  by  absorption  lines  of  OH,  H20, 
CH,  NH3  and  (O  I],  From  the  43-197  pm  continuum  fit 
we  derive  that  the  optical  depth  of  the  thermal  dust 
emission,  xdusti  is  about  unity  at  150  pm.  From  our 
continuum  fit  we  derive  a  Lyman  continuum  production 
rate  of  Q  =  9  x  1054  sec'1  and  our  simple  spherical 
radiative  transfer  solution  indicates  a  cool  central 
source.  The  upper  limits  on  the  FIR  fine-structure  lines 
indicate  a  softer  radiation  field  in  Arp  220  than  in 
starburst  galaxies  such  as  M82  or  in  AGN.  The  low 
L[ciii  is?  /  LFir  ratio  in  Arp  220  cannot  be  easily 
explained  by  dust  obscuration  or  saturation  effects. 
These  results  may  be  explained  by  a  starburst  model 
with  a  low  upper  mass  cutoff  or  by  unusually  high  dust 
absorption  of  UV  photons  within  the  ionized  regions  of 
a  starburst  or  AGN. 

Keywords:  active  galaxies;  infrared  spectroscopy. 

1.  INTRODUCTION 

As  the  closest  (77  Mpc,  H0  =  75  km  s'1)  and  brightest  of 
the  class  of  ultraluminous  infrared  galaxies  (ULIGs), 
Arp  220  (L|R  ~  1.5  xlO12  L0)  is  often  considered  to  be 
the  prototype  of  its  class.  As  is  typical  for  galaxies  in 


this  class,  Arp  220  is  in  the  final  stages  of  the  merging 
of  two  galaxies.  Whether  or  not  Arp  220  contains  a 
luminous  active  galactic  nucleus  that  powers  most  of  its 
infrared  luminosity,  is  a  question  hotly  debated  and 
important  for  our  understanding  of  the  evolution  of 
galaxies  and  of  the  origin  of  active  galactic  nuclei.  We 
try  to  answer  this  question  by  probing  the  hardness  of 
the  radiation  field  using  diagnostic  mid-  and  far-infrared 
fine-structure  lines.  Both  the  detection/non-detection  of 
the  presence  of  an  active  galactic  nucleus  (or  a  young 
starburst)  and  the  derivation  of  the  total  luminosity  of 
these  sources  can  be  highly  sensitive  to  extinction.  For 
example,  if  the  obscuring  dust  is  mixed  with  the 
emitting  gas,  foreground  sources  may  dominate  short 
wavelength  observations  while  the  source  responsible 
for  powering  most  of  the  emission  can  be  effectively 
hidden.  Thus,  to  understand  the  energetics  of  a  galaxy 
as  highly  obscured  as  Arp  220,  far-infrared 
spectroscopy  is  of  critical  importance. 

2.  THE  OBSERVATIONS 

As  part  of  the  LWS  consortium  “Central  Program” 
study  of  infrared-bright  galaxies,  the  observations  of 
Arp  220  and  other  bright  program  galaxies  were 
designed  to  achieve  sufficient  signal-to-noise  to  detect  a 
set  of  seven  diagnostic  far-infrared  fine-structure  lines, 
assuming  line-to-continuum  ratios  similar  to  those  seen 
in  the  nearby  starburst  galaxy  M82  (Lord  et  al.  1996): 
the  [O  III]  52  and  88  pm,  [N  III]  57  pm,  [N  II]  122  pm, 
[O  I]  63  and  146  urn,  and  [C  II]  158  pm  lines.  Over 
13,000  seconds  of  ISO  observing  time  was  devoted  to 
obtaining  a  deep,  full  range  ISO  LWS  43  -  197  pm. 


'  Based  on  observations  with  ISO,  an  ESA  project  with  instruments  funded  by  ESA  Member  States  (especially  the  PI 
countries:  France,  Germany,  the  Netherlands  and  the  United  Kingdom)  with  the  participation  of  ISAS  and  NASA. 
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low  resolution  (R  =  A./AA.  ~  200)  spectrum  of  Arp  220. 
The  LWS  aperture  size  is  ~  70",  encompassing  all  of  the 
far-infrared  flux  of  Arp  220  (Clegg  et  al.  1996; 
Swinyard  et  al.  1996). 

The  data  were  deglitched  and  calibrated  with  the  ISO 
LWS  Pipeline  5.  The  ISO  Spectral  Analysis  Package 
(ISAP)  was  used  to  further  deglitch  and  coadd  the  data. 
The  beginning  and  ending  data  within  each  detector 
scan  were  found  to  be  less  reliable  than  the  central 
regions  of  the  scan.  However,  because  of  the 
redundancy  from  neighboring  detectors,  these  data 
could  be  omitted  still  leaving  sufficient  overlapping  data 
to  be  used  for  scaling  factors  of  <  25%,  our  estimated 
calibration  error,  to  produce  a  continuous  spectrum. 

The  full  LWS  spectrum  of  Arp  220  is  shown  in  Figure 
1 .  Although  the  expected  signal-to-noise  was  achieved, 
with  the  exception  of  weak  [C  II]  emission,  no  fine- 
structure  lines  were  detected  in  emission.  Instead,  the 
spectrum  is  dominated  by  absorption  lines  of  OH,  H20, 
CH,  NH3,  and  [O  I]  and  other  (as  yet)  unidentified  lines. 
For  example,  six  absorption  lines  and  one  emission  line 
from  the  2n1/2,  and  cross-rotational  ladders  of  OH 
and  thirteen  absorption  lines  of  H20,  from  lower  energy 
levels  up  to  300  cm'1,  were  detected.  The  spectrum  is 
markedly  different  from  other  IR  bright,  less  luminous 
starburst  galaxies  in  our  survey,  such  as  M  82,  NGC 
253,  NGC  6240,  and  NGC  3690,  whose  far-infrared 


spectra  are  dominated  by  strong  forbidden  fine-structure 
line  emission,  but  similar  to  that  of  the  ULIG  Mkn  23 1 . 
Here  we  present  the  current  status  of  analysis  of  the 
continuum  and  fine-structure  lines.  Detailed  analysis  of 
the  molecular  lines  will  be  discussed  in  a  future  paper. 


3.  THE  DUST  OPTICAL  DEPTH  AND  THE 
CENTRAL  SOURCE  TEMPERATURE 

To  estimate  the  dust  optical  depth,  xdust(?c),  we  fit  the 
LWS  far-infrared  continuum  spectrum  with  a  function 
of  the  form 

F,<xBx(T)(l-e'TdustW) 

assuming  constant  dust  temperature  and  Tdu,t(%)  = 
(3-oA)a  where  a=  1  -  2,  and  A,0  is  the  wavelength  where 
Tdust  =  L  The  best  fit,  shown  in  Figure  2,  was  obtained 
with  a  =  1.75,  T  =  56  K,  and  \0  =  150  pm.  More 
realistically,  we  fit  the  data  with  a  central  source  +  dust 
shell  radiative  transfer  model  using  a  silicate  (53%)  and 
graphite  (47%)  grain  mixture  (Draine  and  Lee  1984) 

and  allowing  temperature  and  density  (n  cc  r'IJ) 
variation.  We  obtained  a  good  fit  using  two  cool  central 
source  temperatures  (Ti=  5900  K  and  T2=l  15  K),  inner 
and  outer  radii  of  40  and  200  pc  and  dust  temperatures 
of  250  and  50  K,  respectively.  We  show  the  dust 


Figure  1.  The  ISO  LWS  spectrum  of  Arp  220.  The  positions  of  the  strongest  lines  are  indicated. 
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Wavelength  (microns) 

Figure  2.  Best  single  temperature  fit  to  the  Arp  220 
continuum. 


optical  depth  as  a  function  of  wavelength  for  this  model 
in  Figure  3.  Important  for  models  of  the  molecular 
absorption,  is  the  fact  that  the  dust  temperature  at 
"tdustW  ~  2/3  varies  from  TM(180  pm)  ~  220  K  to 
T2/3(20  pm)  ~  60  K.  The  high  values  of  "tdusiM  we 
derive  are  consistent  with  the  high  column  densities  of 
gas  derived  from  millimeter  observations  of  CO 
(Scoville  et  al.  1997).  Any  source  embedded  deep 
within  such  a  dust  shell  may  well  be  hidden  from  view 
even  at  mid-infrared  wavelengths. 


4.  THE  LYMAN  CONTINUUM  RATE  AND  THE 
HARDNESS  OF  THE  RADIATION  FIELD 

To  estimate  the  Lyman  continuum  production  rate,  we 
extrapolated  our  thermal  dust  continuum  fit  to  2.7  mm. 
Subtracting  thermal  dust  and  non-thermai  radio  (Norris 
1988)  extrapolated  components  from  the  2.7  mm 
continuum  flux  (Scoville  et  al.  1991),  we  derive  a 
residual  flux  of  12.6  mJy.  Under  the  assumption  that 
the  residual  is  due  to  optically  thin  free-free  emission, 
we  estimate  a  Lyman  continuum  production  rate  of  Q  = 
8.9  x  1054  s'1.  Comparing  this  to  the  Bra  line  flux 
measured  with  the  ISO  SWS  (Sturm  et  al.  1996)  would 
imply  an  extinction  to  the  starburst  of  Av  =  24 
(assuming  a=1.75)  and  considerably  less  dust  along  the 
line  of  sight  than  implied  by  the  optical  depth  we  derive 
for  the  entire  column  of  dust  (see  Figure  3).  This  could 
be  because  the  starbursts  (at  least  those  we  detect  in 
Bra)  really  do  lie  toward  the  outer  regions  of  the  dusty 
nuclear  regions  or  because  of  the  effects  of  scattering  at 
Bra,  which  we  have  not  taken  into  account. 

To  probe  the  hardness  of  the  radiation  field  of  the 
ionizing  source,  our  upper  limit  on  the  [O  III]  52  pm 


Wavelength  (|im) 

Figure  3.  Dust  optical  depth  as  a  function  of  wave¬ 
length  in  Arp  220  based  on  a  spherical  dust  shell  + 
central  blackbody  model  (see  text). 

(Ep=35  eV)  line  (<7  x  10'20  W  cm'2,  3a)  can  be 
compared  with  the  observed  flux  in  the  [S  III]  34  pm 
(Ep=23  eV)  line  (Sturm  et  al.  1996),  assuming  solar 
abundances.  These  lines  have  nearly  identical  critical 
densities  (ncrjt  -  5  xlO3  cm'3)  and  this  line  flux  ratio  is 
relatively  insensitive  to  the  extinction.  While  in  M82, 
Lord  et  al.  (1996)  observe  [S  III]34  pm  /  [O  III]52  pm  = 
1.35,  the  ISO  observations  of  Arp  220  yield  an  observed 
ratio  of  >  2  (3a).  Higher  extinction  in  Arp  220  makes 
this  difference  more  pronounced.  These  results  would 
appear  to  rule  out  an  AGN  as  the  predominant  energy 
source  for  the  FIR  luminosity  in  Arp  220  or,  if  a 
starburst  powers  the  FIR  emission,  they  imply  that  it 
must  be  older  or  have  a  lower  upper  mass  cutoff  than  in 
M82  (unless  the  central  source  is  hidden  even  at  52  pm 
or  if  dust  in  the  ionized  region  softens  the  spectrum). 


6.  STARBURST  MODELS 

To  derive  the  parameters  of  the  starburst,  an 
evolutionary  synthesis  code  (Leitherer  &  Heckman 
1995)  was  used  to  compute  the  starburst  spectral  energy 
distribution  for  input  to  CLOUDY  (Ferland  1993).  For 
an  instantaneous  burst  with  solar  abundances  and  a 
Salpeter  IMF  (a=-2.35)  for  masses  1  -  100  M0,  we  used 
the  model  [SIII]34  pm  /[OIII]52  pm  line  ratio  as  a 
function  of  age  to  date  the  starburst  using  a  radius  of 
235  pc  (Scoville  et  al.  1997).  Using  the  “Orion  HII 
region  model”  in  CLOUDY,  and  assuming  that  the 
starburst  powers  the  FIR,  a  starburst  age  of  >8x10 
years  and  nc  =  500  cm'3  are  consistent  with  the  observed 
line  ratio  and  line  fluxes,  but  the  Lyman  continuum  flux 
is  a  factor  of  3  lower  than  our  derived  value.  Models 
with  unusually  large  dust-to-gas  ratios  may  allow  a 
better  fit  to  the  observations  (Satyapal  et  al.  1998). 
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[Oil]  15fl  /  FIR 


log  n  (cm'*) 

Figure  4a.  Contour  plot  of  L[c  miss  /LFir  as  a  function 
of  G0  and  n(cm"3)  with  Orion-type  abundances  based 
on  the  models  of  Kaufman,  Wolfire  &  Luhman  (1998). 
For  Arp  220,  this  ratio  is  10'4. 


7.  PHOTODISSOCIATION  REGION  MODELS 

As  in  other  ULIGs  (Luhman  et  al.  1998),  L(c  nnss  /LFIR 
is  lower  by  a  factor  -  10  in  Arp  220  relative  to  typical 
starburst  galaxies.  Since  the  [O  I]  63  pm  line  is  self- 
absorbed,  we  use  our  upper  limit  on  the  [O  I]  145  pm 
line  emission  to  compare  with  PDR  models  to  try  to 
deduce  the  average  incident  far-ultraviolet  radiation 
field  flux  G„  (in  units  of  the  local  interstellar  radiation 
field)  and  PDR  gas  density  n(cm"3).  Figure  4  displays 
L[c  ii] i  5s  4m  /Lfir  and  F(c  11115841m  /F(o  11111454m  as  a  function 
of  n  and  Ga  as  predicted  by  Kaufman,  Wolfire  & 
Luhman  (1998)  in  their  recent  models  including  the 
effects  of  PAH  heating.  These  models  also  include  the 
effects  of  saturation  of  the  [C  II]  line  at  high  n  and/or  Ga 
(Luhman  et  al.  1998).  Based  on  our  [C  II]  158  pm  and 
[O  I]  145  pm  line  flux  measurements  of  8.7  ±  0.04  and 
<  1  (3  a)  xlO"20  W  cm"2  these  ratios  are  10"4  and  a  10 
respectively  in  Arp  220.  Even  if  the  [O  I]  145  pm  line 
flux  is  just  below  our  detection  limit,  from  which  we 
would  derive  G0  ~  20  and  nPDR  -  3  xlO5  cm"3,  the 
derived  value  of  G„  is  unusually  low  for  both  galactic 
star  forming  regions  and  starburst  galaxies.  It  does  not 
appear  that  the  low  [C  II]  line  flux  is  due  to  extinction 
by  dust  at  158  pm,  given  our  derived  Av  to  the  starburst 
and  even  the  total  dust  extinction  through  the  cloud  (see 
§3).  Self-absorption  in  the  [C  II]  line  is  a  possible 
explanation  which  we  are  exploring  with  follow-up 
LWS  Fabry-Perot  observations  (R  =  X  /  AX  ~  104). 
Alternatively,  the  low  [C  II]  line  flux  could  indicate  a 
stellar  population  deficient  in  massive  main  sequence 
stars,  due  either  to  age  or  to  a  low  upper  mass  cutoff 
(see  Satyapal  et  al.  1998). 


[CIIJ  158/(011146 


log  n  (Ctrl'*) 


Figure  4b.  Contour  plot  of  F|c  IIU58tml  /F[0  11111454m 
as  in  Figure  4a.  For  Arp  220,  this  ratio  is  >  10. 


8.  SUMMARY 

We  have  analysed  the  shape  of  the  far-infrared 
continuum  and  the  weakness  of  the  far-infrared  fine- 
structure  lines  in  Arp  220.  Based  on  a  spherically 
symmetric  model  we  find  that  the  far-infrared  emission 
comes  from  a  region  extending  radially  from  40  -  200 
pc  with  Tdust(150  pm)  -  1.  We  derive  a  Lyman 
continuum  rate  Q  ~  9  x  1034  s"1  and  extinction  to  the 
starburst  of  Av  =  24.  The  upper  limits  on  the  FIR  fine- 
structure  lines  indicate  a  softer  radiation  field  in  Arp 
220  than  in  starburst  galaxies  such  as  M82  or  in  AGN. 
These  results  cannot  be  easily  explained  by  dust 
obscuration,  but  may  be  explained  by  an  older  starburst, 
one  with  a  low  upper  mass  cutoff,  or  by  unusually  high 
dust  absorption  of  UV  photons  within  the  ionized 
regions  of  a  starburst  or  AGN. 

This  work  was  made  possible  by  the  hard  work  and 
dedication  of  the  ISO  team  and  the  LWS  consortium. 
We  gratefully  acknowledge  the  use  of  the  Kaufman, 
Wolfire,  and  Luhman  (1998)  models.  This  work  was 
supported  by  the  Office  of  Naval  Research  and  NASA. 
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Question  from  Pierre  Cox 

I 

0 

1 

i  Question:  In  the  case  of  SgrB2,  LWS  Fabry-Perot  mea¬ 

surements  show  that  the  [C  II]  158  [L m  line  is  seen  in 
i  absorption.  It  is  thus  very  likely  that  this  will  also  occur 

in  Arp  220. 


Answer:  We  will  soon  know  the  answer. 


Question  from  Dieter  Lutz 


Question:  Is  the  relatively  small  G  obtained  from  the 
PDR  analysis  still  consistent  with  the  lower  limit  one 
should  be  able  to  obtain  from  the  luminosity  and  size 
of  the  Arp  220  starburst? 


Answer:  G  obtained  in  this  manner  is  on  the  order  of 
105,  in  contrast  to  the  small  G  obtained  by  comparison 
of  the  [C  II]  158/zm  and  [O  I]  145  fjL m  line  observations 
to  the  models.  We  do  not  think  these  lines  are  greatly 
affected  by  extinction  to  the  PDR  regions,  considering 
the  extinction  we  have  derived  to  the  HII  regions.  Per¬ 
haps  these  line  fluxes  are  affected  by  self-absorption.  We 
have  obtained  follow-up  Fabry-Perot  observations  look¬ 
ing  for  self- absorption  in  the  [C  II]  158  /zm  line  in  Arp 
220.  On  the  other  hand,  playing  devil’s  advocate,  we  note 
that  G  can  be  smaller  than  derived  from  this  approxima¬ 
tion  if  the  PDRs  are  located  around  individual  Stromgren 
spheres,  and  not  around  a  starburst  wind,  or  if  the  radi¬ 
ation  field  is  softer  than  typical  of  O/B  stars,  as  we  have 
inferred  from  our  data. 
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